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Abstract 

We demonstrate that the (s-wave) geometric spectrum of the Efimov energy levels in the unitary 
limit is generated by the radial motion of a primitive periodic orbit (and its harmonics) of the 
corresponding classical system. The action of the primitive orbit depends logarithmically on the 
energy. It is shown to be consistent with an inverse-squared radial potential with a lower cut-off 
radius. The lowest-order WKB quantization, including the Langer correction, is shown to reproduce 
the geometric scaling of the energy spectrum. The (WKB) mean-squared radii of the Efimov states 
scale geometrically like the inverse of their energies. The WKB wavefunctions, regularized near the 
classical turning point by Langer's generalized connection formula, are practically indistinguishable 
from the exact wave functions even for the lowest (n = 0) state, apart from a tiny shift of its zeros 
that remains constant for large n. 

PACS: 67.85.-d, 03.65.Sq, 31.15.xg 
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I. INTRODUCTION 



Two particles that are just shy of binding may develop an infinite number of shallow 
bound states when a third particle is added. This was predicted by Efimov [l| forty years 
back, and has only been recently verified experimentally with an ultra-cold gas of optically 
trapped 133 Cs atoms Q]. Subsequently, Barontini et al. have have found evidence for 
two kinds of Efimov trimers with 41 K and 87 Rb atoms. Efimov considered three identical 
bosons interacting pairwise with an interaction whose range tq is much smaller than the 
interatomic scattering length a. Using the hyperspherical coordinates for the three-body 
problem, he showed that the effective potential in the hyperradial coordinate R between the 
length scales of r and \a\ is of an inverse-square type. In the symmetric L = three-body 
state, this effective interaction is sufficiently attractive to give rise to the Efimov spectrum 
for the trimers. A signature of the Efimov spectrum is its geometric scaling, with the ratio 
of the adjacent energy levels being constant. This was predicted by Efimov |l| and has been 
verified by recent experiments 4[]. 

Although in experiments the Efimov spectrum can only be measured for large but finite 
scattering lengths a, our semiclassical analysis of the geometric spectrum in Sec. II is done 
in the unitary limit (\a\ — > oo). The geometric scaling of the spectrum then holds right 
up to the accumulation point at E — 0, and the effective potential is of inverse- square 
type for all distances R > rn. This limit itself has interesting properties and given rise to 
various theoretical research j5KZ[. Our contribution to this research here is a semiclassical 
description of the unitary Efimov system. 

In the periodic orbit theory (POT) there is an intimate connection between classical 
periodic orbits and the quantum energy spectrum of a system through so-called trace formu- 
lae (cf. {9] for an introduction and applications of the POT). In Sec. II, we derive an exact 
trace formula from the Efimov spectrum and show that the action of the classical periodic 
orbit generating the quantum spectrum depends logarithmically on the energy. The corre- 
sponding average level density leads, via an inverse Abel transform, to an inverse-squared 
radial one-body potential that creates the Efimov spectrum in the limit E — > 0. 

In Sec. II we first re-derive the quantum-mechanical solutions of an inverse-squared 
potential with a lower cut-off radius R c . The exact quantum wave function is a modified 
Bessel function of imaginary order {5]. Next we calculate the first-order WKB wave function 
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(including the Langer correction 10j). We show that its leading term in the classically 
allowed region is identical with the leading term of the exact wave function for high-lying 
states (n ^> 1), responsible for the geometric scaling of the spectrum. In the classically 
forbidden region, it decays exponentially like the exact one. The WKB eigenvalues are 
obtained by quantizing the classical action of the inverse-squared potential, whose leading 
term is the action appearing in the semiclassical trace formula discussed in Sec. [Til We find 
that the WKB spectrum, although shlightly phase shifted with respect to the exact one, 
reproduces the geometric scaling of the quantum-mechanical energies and mean-squared 
radii with the same scaling factor. Remarkably, the WKB wave functions, when regularized 
near the classical turning point by Langer's generalized connection formula jloT ]. are - apart 
from their slightly shifted zeros - practically indistinguishable from the quantum-mechanical 
wave functions even for the lowest state (n = 0). 



II. GEOMETRIC SPECTRUM 

In the unitary limit, three identical bosons in the symmetric s state (L = 0), have the 
energy spectrum given by 

E n = E exp(— 2rcn/so) , n — 0, 1, 2, . . . , oo (1) 

so that E n+ i/E n is constant and independent of n and E . This is called the geometric 
spectrum. E < and s > are constants that depend on the system. For three equal- 
mass bosons, one has 

exp(7r/s ) = 22.694 s = 1.00624. (2) 

The energy E , which corresponds to the lowest quantum state of the system (n = 0), 
introduces a length scale whose origin will be made clear soon. Our objective here is to 
derive an exact semiclassical trace formula for the density of states corresponding to the 
spectrum given by Eq. ([T|). This will enable us to identify the action of a single periodic 
orbit that generates the above spectrum. For an energy spectrum governed by only one 
quantum number, there is a rather simple way of deriving a trace formula. Following {9] 
(Chapter 3.2.2), we write E n = f(n), with degeneracy D(n) = 1. The function f(n) is 
monotonic with a differentiable inverse, f~ l {x) = E(x). The exact density of states, defined 
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by 

g(E) = ZJ(E - E n ) , (3) 
can then be rewritten, using Poisson resummation, as 

g{E) = \F'{E)\Y% =G 8{n-F{E)) = \F'(E)\ [1 + 2 E^ =1 cos(2tt kF(E))}. (4) 

This result, which is exact, can be split into two parts: a Thomas-Fermi (TF) term g(E) = 
\F'(E)\ which gives its average behaviour, and an oscillating term which we denote by Sg(E). 

In the semiclassical POT jg], the oscillating part of the exact density of states of a 
quantum Hamiltonian is expressed as a sum over the periodic orbits of the corresponding 
classical Hamiltonian: 



\s v (E)-a T /- 



(5) 



S 9( E ) = ^XMr,fcCos 
r k=i 

The sum is over primitive periodic orbits (k = 1), denoted by T, and their repetitions 
(harmonics) k > 1. The amplitude Ar,k of a periodic orbit (assumed here to be isolated in 
phase space) depends on its primitive period and on its stability matrix. The phase factor 
<Tr,fc is called the Maslov index of the periodic orbit. Comparing the simple trace formula 
(jl]) with the general form ([5]), we see that it has only one primitive periodic orbit T with 
action S(E) = 27ihF(E), and an amplitude 2F'(E) which is proportional to its period. The 
Maslov index is zero. Since the spectrum is for s states only, the action is that of the radial 
motion in a central potential. For the geometric spectrum, inverting (CQ), we have 

n(E) = HE /E) = F(E) . (6) 

Therefore the action S(E) of the primitive orbit is given by 

S(E) = hs \n(E /E). (7) 

Substituting F(E) in Eq. (j4j), we obtain 

9{E) = 2^\E\ j 1 + 2 E cos ^ln(i?o/^)]|. (8) 

This is an exact trace formula, representing a Fourier decomposition of ((3]). When summed 
over all harmonics, it reproduces the quantized spectrum (Tj[|). 
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The smooth part of the density of states g(E) is given by the first term on the r.h.s. of 
OH]). The total number of Efimov states between the energies — \Eq\ and — \E\ is given by its 
integration over this interval, which yields Eq. (jH]). 

Efimov trimers are formed when the two-body scattering length \a\ ^> r , where the 
latter is of the order of the range of the intermolecular potential. One may then take the 
shallowest state to be E a ~ h 2 /ma 2 , and the deepest state to have energy Eq ~ Y?jmr\. As 
a — > ±00, there is an accumulation of states near zero energy. Substituting these in ([6]), we 
find the total number of Efimov states to be 

s 



N 



7T 



ln(|a|/r ) 



(9) 



The same result was obtained by Efimov following a different route. 

In order to obtain a hint to the kind of one-body potential that can generate a spectrum 
of the form (TjQ), we apply an idea ll| that exploits the properties of the Abel transform 



121 ] . The TF level density obtained above, interpreted as the s-state level density of a 



radially-symmetric potential V(r), can be written [9,ll3| as 



9(E) 



so 



m 1 



r-2 



dr 



(10) 



2txE V 2 Tm J n - V(r) 

where r\ and r2 are the lower and upper turning points with V(r±) = Eq and V{r2) = E, 
respectively, and the potential is assumed to be monotonously increasing from its minimum 
value E to the energy E. Using the substitution dr = y(V) dV, so that y(V) = 1/V'(r), 
we can rewrite (TTOj) as 

E y(V)dV 



9(E) 




(11) 



The r.h.s. above represents an Abel transform [12( of the function y(V) which, by the inverse 
transform, is given by 

g(E ) 



y(v) 



J- 

m 



V ~u^ dE 



V^v^Eq 

Using g'(E) = sq/27tE 2 , the integral above is elementary and leads to 

hs 



y(v) 



ixylrn 



+ 



1 



arctan 



E n 



(12) 



(13) 



-v)VV^e-q (-Vf/ 2 V (-V) 

Assuming that V(r) — > for E — > 0, the second term will be leading. We thus get asymp 
totically 



y(v) -> 



hsQ 

2v / 2^ 



-V nj 



-3/2 



VL(r) 



for E . 



(14) 



(Note that both E and V are always negative.) The potential V as (r) which solves the r.h.s. 
above is given by 

U2 2 

An inverse- squared potential of the type (fT5|) is, indeed, shown in the following section 
to be responsible for the asymptotic quantum-mechanical s-state spectrum of the Efimov 
three-body system in the limit E — > 0. 



III. THE INVERSE-SQUARED POTENTIAL 

The three-body problem, after eliminating the center-of-mass coordinates, contains six 
degrees of freedom. In the hyperspherical formalism, these are described by a hyperradius 



R 



(16) 



and five hyperangular coordinates 6j. In the adiabatic approximation, for fixed R, the 
Schrodinger equation for the angular coordinates is solved to obtain a complete set of adia- 
batic eigenstates and the corresponding eigenvalues e(R). In the unitary limit a — > ±oo, the 
angular variables decouple and one gets an effective inverse- squared potential in the coordi- 
nate R. The uncoupled hyperradial wavefunction in the L = state is ^(R) = R~ 5 ^ 2 u(R). 
The reduced wavefunction u(R) obeys the Schrodinger equation (cf. also a pedagogical re- 



view of the Efimov effect 



13|) 



1/4) 



dR 2 



R 2 



u{R) 



2mE 

h 2 



u(R) 



(17) 



where m is the mass of the atom and sq the constant given in (T5]). In order to regularize the 
inverse-squared potential 

K 0( K) *M±IM, m 



2m R 2 

we introduce a lower cut-off radius R c ~ tq as shown schematically in Fig. 1. 
For future convenience we introduce a dimensionless scaled variable x by 



x = R/R A 



x c — R c /R_ 



(19) 



where R + is the classical turning point at the energy E 



E = V(R+) = -hfs 2 j2mR% , 



(20) 
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FIG. 1: A schematic of the inverse-squared potential is shown. R c is the cut-off radius; R + denotes 
the outer classical turning point which is determined by the energy through E = V(R+). 

(see Sec. IIIIBI below), and the lowest energy E c of the classical particle is given by 

E c = V(R C ) = -h 2 sl/2mR 2 c , (21) 

so that 

x\ = E/E c . (22) 



A. Quantum- mechanical solutions 

We now solve the Schrodinger equation ( jTTl) with the lower boundary condition u(R c ) = 0. 
The second boundary condition comes from the requirement that the wave functions vanish 
at infinity. We introduce the following transformation 

u(R) = VRw(R) (23) 

and obtain the equation for the function w(R) as 

d 2 w 1 dw Sn 2mE . . 

IlP + RdR + W W = ~ W - (M) 

Using the scaled variable x = R/R + = R^ —2mE / h 2 , the above equation reduces to the 

standard form for the modified Bessel function 

d 2 w 1 dw (iso) 2 



dx 2 x dx x 2 
7 



w = w . (25) 



Thus, the energy scales away, which is a unique feature of the inverse- squared potential 
whose R dependence is the same as that of the kinetic energy operator. The solution to the 

5I , is a modified Bessel function with imaginary index iso 

w(x) = K iso (y), y = s x , (26) 



equation (J25J) , given m 



which for real y and So is a real function that vanishes exponentially as y — > 00. The 
eigenstates are found from its zeros as explained below. 

n 

To compute the function K iu (y), we use a power series expansion given in [14J j: 
where the phase 4>v,k is given by 

fc 

4> Ui k = arg T(l + k + iv) = <p ufi + ^ arctan ( - ) . (28) 

s =i s 



o and a convergent series for its calculation [15| are given by 

VjO = arg T(l + iv) = -v 7 + V ( — arctan — — ) , (29) 



s=0 

where 7 = 0.577215664 ... is Euler's constant. Numerically we find 

SoO = -0.30103393... (30) 

As stated in [14J], the function K^(y) has an infinite sequence of non-degenerate zeros y n 
(n = 1, 2, . . . ) with < • • • < y n+ \ < y n < y n -i ■ • • < Hi < v , an d no zeros for y>v. The 
solution to the Schrodinger equation ( fl7l) is therefore given by 

u(x) = C<JyK iS0 (y) , (y = s x) (31) 

where C is a normalisation factor. The eigenspectrum is obtained from the zeros y n via the 
boundary condition K iu (y n ) — (n = 1, 2, . . . ). 

It is an important feature of this system that all eigenfunctions u n {x) are given in terms 
of the universal solution fl3T|) simply by letting the variable y in fl3TT) start at the (n+l)-th 
zero y n -\.\. Fig. 2 shows this universal function in a doubly logarithmic plot, exhibiting the 
first 20 zeros. For x > 1, i.e. ln(x) > 0, we notice the beginning of the exponential tail. 
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FIG. 2: The universal function it(x) in a doubly logarithmic plot. We see the first 20 zeros as 
negative spikes for ln(x) < 0. 



We can now associate the zeros with eigenvalues of the scaled cut-off x c and write the 
eigenfunctions as 



Un{%) — C n yf SqX K isQ (SqX) 



X > (x r 



Vn+ll s 



n 



0,1, 



(32) 



where C n is the normalization constant of the n-th state. Note that n here is the number 
of zeros x n j > (x c ) n (j = 0, 1, . . . , n) of the functions u n (x). We found numerically that the 
C n are practically identical for all n > 2 (see Tab. [J below), which is due to the fact that all 
wave functions are peaked around x — 1 and the regions below the first three zeros, y < y^, 
give only exponentially small contributions to the norms. 

In the region of the maximum of u(x) and the exponential tail for y > 1, we had to 
include contributions to the sum in f )27|) up to k max > 25. For y > 11, we did not, in fact, 
obtain convergence of the k sum. However, we found numerically that the leading term of 
an asymptotic form of Ki SQ (y) for large y, given in [l^] 



n {x) ~ C n exp(-y) , (y>l) 



(33) 



becomes sufficiently accurate for y > 9. 
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For y < 0.3, the terms with k > of the series in ( 12?]) become numerically insignificant. 
Since the largest zero is found at y\ = 0.0653423 [cf. ha(x c )o = —2.73434 in Tab. UJ below], 
this means that all zeros of Ki U {y) are given by the leading term with k = 0, yielding the 
asymptotic solution given in |6[ 

K iv (y) cx y/ysm[v\n{y/2)-(j) ViQ ]. (n > 1) (34) 

Its zeros give the geometrical spectrum 

E n /E c = (x c ) 2 n ~ ^^exp(-27rn/so-27r/s o + 20 SOi o/so) (n > 1) (35) 

which, apart from constants, is the same as that discussed in Sec. [Til m particular, the 
constant ratio E n+ i/E n remains the same. From ( 135]) we get 

ln(x c ) n ~ -nir/so + a , (n > 1) (36) 

where a is, using the actual constants in (j2J) and (I3TJ]) . given by 

«o = -vr/so + <p S(h o/s + hi(2/s ) = -2.73434 . (37) 

In column 1 of Tab. [TJ we give some selected eigenvalues ln(x c ) n obtained numerically 
from the exact solutions above. Although the expression (136]) is mathematically correct 
only asymptotically for large n, we find that the exact numerical eigenvalues agree with the 
r.h.s. of ( 136]) and ( 137]) within five digits - which corresponds to our numerical accuracy - 
even for n = 0. 



n 


ln(x c ) n (QM) 


In(x c )„ (WKB) 


5\n(x c ) n 


1/Cn 





-2.73434 


-2.64717 


0.08717 


0.3649953 


1 


-5.85644 


-5.77053 


0.08591 


0.3651889 


2 


-8.97854 


-8.89263 


0.08591 


0.3651892 


10 


-33.95535 


-33.86943 


0.08592 


0.3651892 


11 


-37.07745 


-36.99154 


0.08591 


0.3651892 


20 


-65.17635 


-65.09044 


0.08591 


0.3651892 


63 


-199.42668 


-199.34078 


0.08590 


0.3651892 


100 


-314.94440 


-314.85850 


0.08590 


0.3651892 



TABLE I: Eigenvalues ln(x c ) n of some selected states, obtained numerically both quantum- 
mechanically (QM) and in the WKB approximation discussed in Sec. IIIIBI Column 3 gives 
their differences, and column 4 shows the inverse normalization constants C n in (|32p . 
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It is quite remarkable that, with the value so — 1 given for the present system, the 
geometrical nature of the Efimov spectrum is preserved with a high numerical accuracy all 
the way down to n = in the inverse- squared potential. For appreciably larger values of s 
this would not be the case. 

We emphasize once more that the geometric nature of the asymptotic Efimov spectrum, 
given by equations ([1]) and ( )35|) . is a direct consequence of the inverse- square behavior of 
the potential (118)) . with its strength determined by the parameter sq. The ground-state 
energy Eq in ([T|), and with it the parameter <fr So ,o m (|35|) . is determined by the particular 
choice of the regularization of the potential at small distances. In the present case, this 
was done by the hard-wall cut-off at R c and the corresponding boundary condition. Other 
forms of regularization would lead to the same geometric scaling as in ([T]) but to different 



ground-state energies E , such as e.g. in 
addition of a repulsive 1/R A potential at short distances. 



161 ] where the potential f fT8|) was regularized by the 



B. WKB solutions 



We now present our results using the WKB approximation. Since the WKB method is 



well known, we need not present it here and refer for details to the text book by Migdal 17| . 



who also discusses explicitly the Langer correction for radially symmetric potentials [10 1. 

In the scaled variable x, the classical orbits are bounded inside the classically allowed 
region x c < x < 1. Applying the Langer correction with the potential ( !T8|) . the effective 
classical potential becomes 

v (R) = -^§, (R>R C ) (38) 

which is identical with the asymptotic potential ( I15p found from the inverse Abel transform 
of the TF level density in Sect. II. The classical momentum then is 



P(R) = \J2mE + h 2 s 2 /R 2 , (39) 
which, in the scaled variable x, becomes 

P(R) = ^ = ~V0^j, Po = ^- (x c <x<l) (40) 
x R+ 



The standard WKB wave function in the classically allowed region has the form 

(x c < x < 1) (41) 



WKB / \ A 

uZ [x) = -p= cos 



j_S in (x) - 
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where Si n (x) is the action integral along the orbit from x to 1: 



S in (x) = / p(x)dx = hs 

J x 

This integral can be found analytically 

S in (x) = hs 



x- 



dx 



x 



(x c < x < 1) 



(42) 



and becomes 



In + \n(l + Vl - £ 2 ) - VT 



(x r < X < 1) 



(43) 



The boundary condition that the wave function (14"T!) vanishes at the lower turning point 
x = x c leads, like in the quantum-mechanical case, to the quantization of the eigenenergies 
as shown below. 

Outside the classically allowed region (x > 1), the exponentially decaying WKB wave 
function has the form 



B 



exp 



y/K(x) 

with the action S out (x) given by 



'^:S ou t(x) 
n 



K {x) = ^ViM, (x > 1) (44) 



Sout(x) = I k(x) dx = hs 



arctan 



Vx^l 



+ Vx 2 - 1 - - 



(x > V 



(45) 



For large x 3> 1, the function (J44J) has - apart from the normalization - the same exponential 
tail as the exact function (1331): 



(46) 



Equating the two asymptotic forms (1331) and (1461) . we can determine the constant B as 



Be n/2 

u ™f B ( x ) ^=ir exp(-s x) . (s x > 1) 



'SqK 



B n = C n ^s h7i/2e- 7r 



/2 



(47) 



so that both wave functions agree exactly in the asymptotic tail region for each state \n). 

Since both WKB functions ( )4T|) and ( 1441) diverge at the upper classical turning point 
x = 1, where p(l) = k(1) = 0, one must regularize these functions. This is done in the 



standard "connection" 
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17J by a linear approximation to the potential V(R) in the 



neighborhood of R = R + and matching the corresponding Airy function solution to the 
above WKB wave functions at some distances on either side of x = 1. The requirement that 
the wavefunction be continuous and continuously different iable at both matching points 
leads to the relation A n = 2B n for the normalization constants and the phase —ir/4 in (j4T|) . 
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We now have to fulfil the lower boundary condition at x = x c : 

= 0, (48) 



.wkb,„ ^ _ ^ CQS 



where Sq(E) is the action of the primitive periodic orbit going from x c to x = 1 and back: 

S (E) =2/ p(z)dz, (49) 



which, using fl43|) and (j22|) becomes 

S (E) = hs hx(E c /E) + 2 In (l + y/l - E/E^j - 2y/l - £/£ c . (50) 

For <C |i? c |, i.e. for the the shallow states (n ^> 1), the second and third terms above 
become negligible and the leading term reproduces the action Q (with E = E c ) obtained 
from the periodic orbit theory applied to the geometric spectrum, as discussed in Sec. [Ill 
Eq. (1481) has the solutions 

L S (E)-^= (n+ 1 )tt, n = 0,l,2,... (51) 



2/i y 4 V 2 
which yields the WKB quantization condition 



SbCE™ ) = J P{R) dR = 2nh (n + 3/4) , n = 0, 1, 2, . . . (52) 

Note that the constant 3/4 in (152]) differs from the usual value 1/2, which one obtains for 
smooth potentials, due to the hard-wall reflection at the lower turning point x c . 

Using only the asymptotically leading logarithmic term in (|50|) yields the geometric spec- 
trum 

e wkb _ Ec exp [_ 27r / So ( n + 4/3)] ; (n > 1) (53) 
which corresponds to 

ln(x c ) n nir/so - 3tt/4s = -nix / s - 2.34158 . (n > 1) (54) 

This is the same as the quantum- mechanical result ( )36l) . apart from a different shift (denoted 
by «o there). In column 2 of Tab. HI we give the WKB eigenvalues obtained from the 
quantization condition ( 1521) . using the full action ( 1501) . in terms of the scaled values ln(x c ) n . 
We see that they come very close to the exact quantum values. In fact, there remains a slight 
shift in ln(x c ) n that becomes constant for n > 1 (cf. the third column in Tab. [I]). A similar 
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(but different) shift between the asymptotic exact and WKB spectra has been obtained for 
an attractive 1/R 2 potential regularized differently 16j (see 19j for its interpretation). 

We now have to regularize the WKB functions (141 p . (1441) near the turning point (x = 1) 
where they diverge. While the standard connection to the Airy solution of the locally 
linearized potential leads to the WKB quantization condition (I52p . as described above, we 
found that it does not yield satisfactory wave functions. The reason is that the asymptotic 
Airy solutions do not come sufficiently close to the WKB solutions on either side of the 
turning point. We were, however, successful when u sing Langer's generalized "connection 
formula" for one isolated classical turning point [Ref. 10], Eqs. (11a) and (lib) with n = 0]. 
Expressable in terms of a single Airy function Ai(£), it reads 



u 



r( X ) = D n v^jTiii Ai(o , x > (x e ), 



(55) 



where D n is a normalization constant, and 

2/3 



2h Sout{ - x) 



2h Sm{x) 



2/3 



Six) 



S(x) 



s, 



out 



X 



hp(x) 



for x > 1 , f > , 
for x < 1 , f < . 



(56) 
(57) 



The superscript "uni" in (1551) stands for "uniform" , because it turns out that fl55|) is a global 
uniform approximation that can be used not only in the vicinity of the classical turning 
point, but throughout the whole domain (x c ) n < x < oo. By construction [10|j, it yields the 
asymptotic WKB solutions (j4ip . (1441) sufficiently far from the turning point x — 1. 

The normalized uniform WKB wavefunction (|55|) for n = 5 is shown in Fig. 3, along with 
the exact one, in a doubly logarithmic plot. Apart from the slightly shifted zeros, there 
is a surprisingly good agreement. Fig. 4 shows the ground-state wave function for n = 0, 
both exactly (solid line), the "raw" WBK approximation (dotted line), and in the uniform 
approximation ( 1551) (dashed line). Even for this lowest state, the quantum-mechanical and 
the uniform WKB solutions are practically indistinguishable. 



Note added after publication of this paper: 

The above global uniform approximation is also discusse in: 
C. M. Bender and S. A. Orszag: Advanced Mathematical Methods for Scientists and Engi- 
neers (Springer- Verlag New York 1991), Ch. 10.4, pp. 510 ff. 
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X 

FIG. 4: The ground state uq(x) in a linear plot. Solid line: exact normalized quantum-mechanical 
solution (|32p . Dotted line: WKB approximation (|4ip . (|44p . normalized to the exact solution in the 
tail region via equation (|47p . Dashed line: normalized uniform WKB approximation (|55p . 
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C. Mean squared radii 



The size of our system in the n-th state is given by the mean squared hyperradius defined 



as 



R 2 ) 





R 2 






n 





R\{En) f 



dxx 2 \u n {x)y 



(58) 



where Eq. (fT9|) and the normalized wavefunctions u n (x) have been used. Since R 2 + scales 
like 1/E n , see (1201) . the ratio of mean squared radii of two successive states is E n /E n+ i, if 
the expectation values (n\x 2 \n) remain independent of n. 

That this is, indeed, the case for all n > 1, is demonstrated in Tab. [Til both for the 
quantum-mechanical and the uniform WKB solutions. Thus the mean squared radius in- 



creases by the same geometric scaling factor as the energy decreases, as discussed in jg], also 
in the WKB approximation. 



n 


{n\x 2 \n) (QM) 


(n|x 2 |n) (WKB) 





1.3265 


1.3408 


1 


1.3251 


1.3392 


2 


1.3251 


1.3392 


5 


1.3251 


1.3392 


10 


1.3251 


1.3392 



TABLE II: Mean squared radii (n\x 2 \n) of some states, obtained quantum- 
mechanically (QM) and in the uniform WKB approximation. 
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IV. SUMMARY 



To summarize, we have shown that the (s-wave) geometric spectrum of the Efimov energy 
levels in the unitary limit is semiclassically generated by a single periodic orbit whose action 
depends logarithmically on the energy. The smooth part of the s-state level density, obtained 
by the periodic orbit theory, is consistent with an attractive inverse-squared radial one-body 
potential. We have re-derived the quantum spectrum of the inverse- squared potential with 
a lower cut-off and shown, for sq — 1, that it reproduces the geometric Efimov spectrum 
not only for shallow states, but yields the same constant ratio E n+1 /E n all the way down 
to the ground state with n = 0. We have given an analytical expression for the zeros of 
the eigenfunctions. The WKB quantization of the classical system (including the Langer 
correction) yields the same spectrum, although slightly phase shifted, which preserves the 
same constant ratio E n+ \/E n down to n — 0. The action of the classical system has as its 
leading term precisely the action obtained from the periodic orbit theory. The WKB wave 
functions, when regularized around the classical turning point using Langer's generalized 
connection formula, reproduce the exact ones surprisingly well, apart from the slightly shifted 
zeros - even for the ground state (n = 0). Both the quantum and the WKB solutions 
reproduce the geometric scaling of the mean squared radii of the Efimov states which is 
inverse to that of their eigenenergies. 

As mentioned in the introduction, the Efimov spectrum is experimentally not measured 
in the unitary limit, but for finite and large scattering lengths \a\ 3> tq. Strictly speak- 
ing, therefore, the 1/R 2 nature of the effective potential is only guaranteed in the range 
R c < R < \a\. Then, the geometric scaling of the spectrum would only hold within the 
corresponding energy range, the approximate number of Efimov bound states being given 
by Eq. 09]). Also, the derivation of the inverse-squared potential given in Sec. II and the 
analytical exact wavefunctions given in Sec. Ill A would no longer hold. However, once we 
assume an inverse-squared potential in the range R c < R < \a\, our WKB calculations still 
go through provided we restrict the outer turning point to R + < \a\. 

R.K.B. and M.B. are grateful to the IMSc, Chennai, for its hospitality, excellent working 
conditions and financial support. We are also grateful to an anonymous referee for valuable 
suggestions. M.B. acknowledges enlightening discussions with H. Friedrich, C. Eltschka and 
D. L. Sprung. 
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